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TESTING  OF  HYPOTHESES  FOR  DISTRIBUTIONS  IN 
ACCELERATED  LIFE  TESTING 
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Moshe  Shaked 
Nozer  D.  Singpurwalla 


1.  INTRODUCTION 

It  is  a  common  practice  to  subject  long-life  items  to  larger 
than  usual  stresses  so  that  failure  data  can  be  obtained  in  a  short 
amount  of  test  time.  Such  tests  are  called  accelerated  life  tests,  and 
the  goal  is  to  infer  the  mean  time  to  failure,  or  more  generally,  the 
life  distribution  of  the  items  at  the  usual  stress  level  using  failure 
data  from  accelerated  tests. 

A  parametric  approach  to  this  problem  [3,5,10]  involves  making  two 
assumptions:  first,  the  life  distribution  at  each  stress  level  is 
assumed  to  come  from  the  same  parametric  family  which  is  specified; 
second,  a  relationship  called  the  "time  transformation  function"  is 
assumed  among  the  parameters  of  the  various  distributions.  As  an  ex¬ 
ample  of  the  former,  it  is  typical  to  assume  that  the  life  distributions 
at  the  various  stress  levels  are  either  exponential,  Weibull,  lognormal, 
or  an  extreme  value.  The  time  transformation  functions  which  are  com¬ 
monly  assumed  are  the  inverse  power  law,  the  Arrhenius  law,  the  Eyring 
law,  or  generalizations  of  these.  Whereas  the  above  assumptions  may  be 
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reasonable  in  many  situations,  the  possibility  does  remain  that  under 
certain  circumstances  an  analyst  may  be  hesitant  or  unwilling  to  enter¬ 
tain  any  one,  or  both,  of  these  assumptions. 

In  this  paper,  we  drop  the  requirement  that  the  common  parametric 
family  of  life  distributions  under  all  stresses  be  specified,  but  retain 
the  requirement  that  the  time  transformation  function  be  specified.  The 
time  transformation  function  that  we  consider  is  a  version  of  the  in¬ 
verse  power  law.  By  choosing  this  form,  we  do  not  sacrifice  any  lack  of 
generality,  since  the  procedures  mentioned  here  can  easily  be  extended 
to  consider  the  other  time  transformation  functions  as  well.  The  in¬ 
verse  power  law  being  commonly  used  is  considered  here  for  purposes  of 
illustration.  We  indicate  how  the  life  distribution  and  the  mean  time 
to  failure  at  use  conditions  stress  can  be  consistently  estimated  and 
how  hypotheses  about  the  common  parametric  family  of  life  distributions 
can  be  tested.  Thus  we  are  able  to  test  the  validity  of  the  distribu¬ 
tional  assumption  traditionally  made  in  accelerated  life  testing,  if  the 
form  of  the  time  transformation  function  can  be  specified.  In  another 
paper  [6],  both  the  assumptions  mentioned  before  have  been  dropped  and 
the  problem  has  been  treated  completely  nonparametrically .  However,  the 
statistical  precision  of  the  results  obtained  there  is  less  than  that 
which  could  be  achieved  by  making  one  or  both  of  the  aforementioned 
assumptions . 

We  should  remark  that  the  dual  problems  of  estimation  and  testing 
for  the  goodness  of  fit  at  use  conditions  stress  have  already  been  con¬ 
sidered  in  two  separate  papers  by  Shaked,  Zimmer,  and  Ball  [9],  and  by 
Sethuraman  and  Singpurwalla  [8],  These  papers  emphasize  technical  de¬ 
tails  pertaining  to  a  formal  justification  of  the  methodology,  the  de¬ 
velopment  of  the  appropriate  formulae,  and  the  performance  of  the  pro¬ 
cedures.  Thus,  a  practitioner  interested  in  using  the  methodologies  may 
find  it  difficult  to  extract  that  information  which  is  pertinent  to  his 
goals.  Furthermore,  the  results  of  the  above  two  papers  complement  each 
other  and  can  therefore  be  combined  to  constitute  a  useful  package  for  a 
nonparametric  approach  to  inference  from  accelerated  life  tests. 

-  2  - 
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Our  goal  In  writing  this  paper  is  to  combine  the  techniques  of 
[8]  and  [9]  in  a  manner  than  is  comprehensive  and  user-oriented  so  that 
this  nonparametric  approach  to  acceierated  life  testing  is  accessible  to 
the  reliability  analyst.  We  demonstrate  the  utility  of  this  approach  by 
applying  it  to  some  data  from  a  realistic  situation. 


2.  A  MODEL  FOR  ACCELERATED  LIFE  TESTING 

Let  V, ,V. . V,  and  V.  denote  the  accelerated  and  the  use 

12  k  0 

conditions  stresses,  respectively,  and  let  F^,...,F^  and  F^  denote 

the  corresponding  cumulative  distribution  functions.  We  shall  assume 
that  F  , ....F^  belong  to  a  common  but  unknown  family  of  distributions, 

and  that  for  some  distribution  F  ,  which  also  belongs  to  this  family, 

F±(t)  =  F(Av“t)  ,  i=0,...,k  ,  (2.1) 

where  A  >  0  and  a  >  0  are  unknown  constants. 

This  assumption  means  that  a  change  of  stress  does  not  change  the 
shape  of  the  life  distribution,  but  only  changes  its  scale.  This  model 
is  a  generalization  of  the  familiar  power  law  model  which  is  tradition¬ 
ally  applied  to  the  scale  parameter  of  the  failure  distribution  [4,  5, 
11]. 


Our  goal  is  to  estimate  F^  ,  the  life  distribution  under  use 

conditions  stress,  and  to  perform  some  goodness  of  fit  tests  for  hypoth¬ 
eses  about  the  general  family  of  distributions  to  which  F  ,  and 
V  ••• » *  belong. 


3.  ESTIMATION  OF  FQ 


Typically,  data 
of  observations  T  ^  , 


from  an  accelerated  life  test  consists  of  the  set 
£=l,...,n^  ,  i=l,...,k  ,  where  T  ^  is  the  time 
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to  failure  of  the  Ith  item  in  a  sample  of  size  n^  that  is  run  under 
the  constant  application  of  stress  ,  i=l,...,k  .  We  assume  that 

,  k  ,  and  n^,...,n^  are  all  fixed  in  advance.  If  data  at 
the  nonaccelerated  stress  level  are  also  available,  then  the  proce¬ 

dure  to  be  described  can  still  be  used,  but  by  augmenting  a  0  to  the 
range  of  indices  i  and  j  below. 


3.1  Background 


Denote  the  scale  factor  between  F.  and  F.  by  0..  ,  where 

i  J  ij 

6..  =  AV°7AVa  =  (V./V.)a  ,  i*j  .  (3.1) 

13  J  i  J  i 

The  first  step  in  the  procedure  for  estimating  F^  is  the  estimation  of 
a  .  A  method  for  estimating  a  is  suggested  by  (3.1),  since 

a  =  log(0  _)/(logV.  -  logVi)  . 


-1  r 

Let  T.  =  n.  I  T ,  i=l,...,k  ,  be  the  sample  means  of  the 
1  1  1=1  1)6 

T  's  ,  S,=l,...,n,  .  Then,  an  estimator  of  0..  is 
i£  ’  i  ij 

8d  ■  Vfj  ■  • 

from  which  we  can  obtain  an  estimator  of  a  as 


n . 

l 


c^.  =  log(O^)/ (logVj  -  logV^)  ,  i^j  . 

The  ^k(k-l)  estimators  of  a  can  be  used  to  form  a  weighted  average 
of  the  a^'s  >  which  gives  us  the  final  estimate  of  a  ,  as 


u 


l  l  (log(V  /V  )J (log(T  /T  ) J 
i=lj=i+l  »  3  '\  J  » 


k  k  i  », 

l  l  log(V  /V  )  ' 
i=l  j=i+l  »  3  1  / 
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If  a  were  known,  we  could,  in  the  light  of  relationship  (2.1), 
transform  (rescale)  the  times  to  failure  from  any  stress  level  to  corres¬ 
pond  to  the  times  to  failure  from  any  other  stress  level.  Thus,  observ¬ 
ing  ,  £=l,...,n^  under  stress  is,  under  our  model,  equivalent 

Cl 

to  observing  (V./V.)  T.„  ,  £=l,...,n.  under  stress  V.  .  Since  a  is 
1  j  xx-  1  j 

not  known,  we  shall  use  its  estimate  a  to  rescale  our  variables.  This 

is  the  key  notion  which  enables  us  to  obtain  an  estimate  of  . 


't.J 

.•’n 

!  I 


3.2  Estimation  of  and  the  Mean  Time 


to  Failure  under  V 


0 


a  vk 

Using  a  given  by  (3.2),  we  define  the  N  =  n^  rescaled 

variables 


U 


'  k)  T±l  ’ 


1, . . .  ,n^  ,  l  I,..., k  , 


(3.3) 


and  estimate  FQ  by  the  empirical  distribution  function  of  the  T^'s  . 
Thus,  Fq  ,  an  estimator  of  the  failure  distribution  at  use  conditions 
stress  Vq  ,  is  given  by 

Ffl(t)  =  ^  (number  of  observations  among  the  T^'s  <  t}  , 


for  0  <  t  <  °°  . 

An  estimator  of  the  mean  time  to  failure  at  use  conditions 
stress,  say  p  ,  is  given  by 


(3.4) 


.  n . 

k  l 


|  l  1  it 

N  i=l  £-1  3 


k  nlv. Is 

i-i  *  Iv  T‘  ' 


(3.5) 


The  estimators  F^  and  p  are  s-consistent  [91.  A  method  for  obtaining 

•A 

the  approximate  s-uniform  confidence  limits  for  F^  using  F^  is  dis¬ 
cussed  in  Section  5. 


fc 
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4.  GOODNESS  OF  FIT  TESTING  OF  HYPOTHESES  FOR  FQ 

Having  estimated  F^  by  F^  given  by  (3.4),  our  next  goal  is  to 

see  if  Fg  is  a  member  of  a  specified  family  of  distribution  functions. 

This  type  of  information  may  be  of  interest  in  its  own  right,  or  it  can 
be  used  to  obtain  s-uniform  confidence  bounds  for  F^  .  Testing  to  see 

if  Fg  is  a  member  of  a  specified  family  is  equivalent  to  testing  wheth¬ 
er  F^,...,F^  are  also  members  of  the  same  family.  For  reasons  given  in 

[8],  it  is  convenient  for  us  to  consider  the  random  variables  X.„  = 

lJt 

logTi^  ,  £=l,...,ni  ,  i=l,...,k  ,  and  test  for  the  underlying  distribu¬ 
tion  of  the  X^'s  >  say  H  ,  i=l,...,k  .  For  example,  if  we  wish  to 
test  the  hypothesis  that  the  F^'s  are  members  of  the  family  of  Weibull 
distributions,  then  this  is  equivalent  to  testing  whether  the  H^'s  are 

members  of  the  family  of  extreme  value  distributions  f 3 ] .  Similarly,  if 
we  wish  to  test  the  hypothesis  that  the  F^'s  are  members  of  the  lognor¬ 
mal  family  of  distributions,  then  this  is  equivalent  to  testing  whether 
the  H.'s  are  normal  [3]. 

l 

There  are  two  limitations  to  our  testing  procedure.  The  first 
one  is  that  we  can  only  work  with  complete  (uncensored)  samples,  and  the 
second  one  is  that  the  H^'s  can  only  be  tested  for  being  members  of  the 

s-location-scale  family  of  distributions.  The  commonly  used  members  of 
the  s-location-scale  family  are  the  normal,  the  exponential,  and  the 
extreme-value  distributions.  Since  H^  is  an  exponential  distribution 

whenever  F.  is  a  Pareto  distribution,  and  since  a  Pareto  distribution 

l 

is  rarely  used  to  describe  life  lengths,  the  desire  for  testing  whether 
H^  is  an  exponential  is  not  very  common. 

The  procedure  for  testing  depends  upon  the  particular  family  of 
distributions  that  is  being  considered;  this  will  be  made  clear  in  the 

-  6  - 
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following  section.  Furthermore,  the  procedure  is  valid  when  each  of  the 
n^'s  is  large, 

4.1  Testing  Whether  is  a  Weibull 


Suppose  that  we  wish  to  test  the  hypothesis  that  the  common  under¬ 
lying  family  of  distributions  is  a  Weibull  distribution.  This  is  equiva¬ 
lent  to  testing  whether  the  H.'s  are  extreme  value  distributions 

1 


H.(x) 


i=l, . . .  ,k 


where  and  6^  are  the  s-location  and  the  s-scale  parameters,  re¬ 

spectively.  The  test  procedure  entails  the  following  steps. 

(1)  Using  ,  £=l,...,n  ,  we  first  obtain  the  s-maximum 

likelihood  estimates  of  fit.  and  3.  ,  say  a.  and  6.  , 

l  l  i  l 

respectively,  using  the  above  specified  .  Methods  for 

obtaining  these  estimates  are  given  in  [3];  they  can  also 
be  obtained  using  the  computer  program  discussed  in  Sec¬ 
tion  4.3. 

(2)  For  0  <  t  <  1  ,  we  compute 


W.(t) 

where 


n . 
i 


L  l  I  HI 


ni  «,=] 


<  t 


and  l{E}  denotes  the  indicator  of  the  event  E  .  Verify 
that  W^(t)  is  merely  the  empirical  distribution  function 


of  the 


Hi((xi{-Si)/31)  ,  £=!,..., n.  . 
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(3)  We  then  compute,  for  0  < t  <  1  , 

V.(t)  =  /Z  jW.(t)-t), 

and  repeat  the  Steps  (1),  (2),  and  (3)  for  all  values  of 
i  ,  I- 1 ,  •  •  • ,  k  « 

(4)  Our  test  statistic  D  is  the  s-Kolmogorov-Smirnov  statis¬ 
tic  obtained  by  pooling  the  V^tJ's  in  a  manner  shown 
below.  Let 


A  =  n.  /  7  n.  ,  and  let 

i=l  1 


Z(t)  =  }'  A.V.(t)  ,  0  <  t  <  1  . 

i=l  1  1 


Then  D  =  max(D  )  ,  where 

i 

X 

1 


D+  =  max  |  y  A2^  Z(t)  ,  and 
0<t<l\i=l 


D~  =  min  (  y  A2)  Z(t)  . 
0<t<l\i=l  Xi 


(5)  We  will  accept  the  hypothesis  that  the  common  underlying 

family  of  distributions  is  a  Weibull  (i.e.,  the  H^'s  are 

extreme  value  distributions)  at  a  desired  level  of  signif¬ 
icance,  say  a  ,  if  the  value  D  is  less  than  the  critical 
values  given  below.  These  values  have  been  abstracted 
f  rom  [21. 


Level  of  Significance,  <x 

.10 

.05 

.025 

Critical  Values  of  D,  D* 

.800 

.870 

.940 
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4.2  Testing  Whether  is  Lognormal 


Testing  to  see  if  the  common  underlying  family  of  distributions 
is  a  lognormal  is  equivalent  to  testing  whether  the  H/s  are  normal 
distributions. 


"i  111  u-a  • 

‘  “H' H"®!1 


i=l,...,k  ,  where  the  and  B^  are  the  s-location  (mean)  and  s- 

scale  (standard  deviation)  parameters,  respectively.  The  test  proce¬ 
dure  entails  the  following  steps. 


(1)  Using  ,  £=!,..., n^  ,  we  obtain 


^i  =  E  Xiff  1  ni  and  =  l  (X  p -at ;  ) 2  /  ("i"1)  . 

1=1  1=1 

the  s-maximum  likelihood  estimators  of  a.  and  B.  , 

l  t 

respectively. 


(2)  We  repeat  steps  (2),  (3),  and  (4)  of  Section  4.1,  except 
that  now 


i  rai 


H.ppA  -  fiZ-±-eJ-i 

\  B.  /  S2v&.  1  2 


S,=l,2, . . .  ,n^  ,  and  to  test  for  the  significance  of  D  ,  we 
use  the  following  critical  values  abstracted  from  [7J. 


Level  of  Significance,  a 

.10 

— 

.05 

.025 

.01 

Critical  Values  of  D,  D* 

.752 

•835j 

.910 

.991 
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4.3  Some  Comments  on  the  Test  Procedure 

The  power  of  the  above  tests  against  various  alternatives  is  dis¬ 
cussed  in  [8],  The  main  point  here  it  that  it  is  difficult  to  distin¬ 
guish  between  the  Weibull  and  the  lognormal  alternatives,  so  that  the 
final  choice  of  the  family  of  distributions  to  which  belongs  must 

be  based  on  criteria  other  than  a  formal  goodness  of  fit  test.  For  ex¬ 
ample,  if  the  hazard  rate  is  known  to  be  increasing  in  time,  then  the 
distribution  cannot  be  lognormal  and  the  Weibull  assumption  should  be 
chosen. 


To  aid  in  the  computation  of  f)  ,  and  also  the  s-maximum  likeli- 
hood  estimates  and  for  the  various  hypotheses  that  we  wish  to 

entertain,  we  have  written  a  computer  program,  described  in  [1].  This 
program  can  be  made  available  to  the  interested  reader. 

5.  UNIFORM  CONFIDENCE  BOUNDS  FOR  Fq 

The  approximate  s-uniform  (l-a)%  confidence  bounds  for  FQ  can 
be  obtained  once  a  test  of  hypothesis  results  in  the  acceptance  of  the 

A 

hypothesis.  The  approximate  confidence  bounds  will  be  based  on  F^ 

given  by  (3.4)  and  the  critical  values  D*  given  in  Section  4.1  and 

4.2.  (The  confidence  bounds  are  approximate  because  of  the  finiteness 

/\ 

of  the  sample  sizes,  and  because  they  are  based  on  and  not  on  a 

complicated  transformation  of  (Z(t),  t >  0}  of  Section  4.)  These 
bounds  can  also  be  used  to  obtain  an  approximate  (l-a)%  s-confidence 
bound  on  u  ,  the  median  of  the  time  to  failure  distribution  under  VR  . 

Suppose  that  the  procedure  of  Section  4.1  results  in  our  accept¬ 
ing  the  hypothesis  that  the  are  extreme  value  distributions.  This 

means  that  FQ  is  a  Weibull  distribution,  so  that  the  approximate  uni¬ 
form  (l-ot)%  confidence  limits  for  F^  are  given  by 
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FQ(t)  ±  D*/^J  , 

where  D*  is  the  critical  value  of  D  corresponding  to  a  ,  and  N  = 

n^  .  For  example,  if  we  choose  a  =  .05  ,  then  the  approximate  95% 
uniform  upper  and  lower  confidence  limits  for  F^  are 

g(t)  =  FQ(t)  +  .870/v^  , 

and 

h(t)  =  FQ(t)  -  . 870/^N  , 

respectively.  In  computing  the  above,  we  should  make  sure  to  confine 
both  g(t)  and  h(t)  to  be  between  0  and  1. 

In  a  similar  manner,  if  the  procedure  of  Section  4.2  results  in 
our  acceptance  of  the  hypothesis  that  F^  is  a  lognormal  distribution, 

then  the  95%  uniform  confidence  bounds  for  Fg(t)  will  be  F^Ct)  - 

.835/^fr  . 

The  uniform  confidence  bounds  g(t)  and  h(t)  can  be  used  to 
obtain  an  approximate  s-confidence  bound  for  jj  by  taking  their  in¬ 
verses.  Thus,  the  approximate  s-confidence  bounds  for  u  are  g  \.5) 
and  h  *(.5)  ,  where  g  ^(.5)  <  h  *(.5)  . 

6.  EXAMPLE 

We  apply  the  methodology  of  this  paper  to  some  accelerated  life 
test  data  given  in  [4]  and  [5].  These  data  represent  the  times  to 
breakdown  of  an  insulating  fluid  subjected  to  seven  elevated  voltage 
levels.  Since  the  methodology  of  Section  4  is  valid  for  large  sample 
sizes,  and  since  three  of  the  stress  levels  of  [5]  contain  too  few 
failures,  we  delete  them  from  our  consideration,  and  consider  an  ab¬ 
stracted  version  of  the  data.  These  data  are  given  in  Table  l,  and  have 
also  been  considered  in  [6].  The  use  conditions  stress  corresponds  to 
Vq  «  20  kilovolts. 
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TABLE  1 


TIMES  TO  BREAKDOWN  OF  AN  INSULATING  FLUID 
(IN  MINUTES)  UNDER  VARIOUS  VALUES  OF  THE 
STRESS 


36  kV 

34  kV 

32  kV 

30  kV 

.35 

.19 

.27 

7.74 

.59 

.78 

.40 

17.05 

.96 

.96 

.69 

20.46 

.99 

1.31 

.79 

21.02 

1.69 

2.78 

2.75 

22.66 

1.97 

3.16 

3.91 

43.40 

2.07 

4.15 

9.88 

47.30 

2.58 

4.67 

13.95 

139.07 

2.71 

4.85 

15.93 

144.12 

2.90 

6.50 

27.80 

175.88 

3.67 

3.99 

7.35 

8.01 

53.24 

82.85 

194.90 

5.35 

8.27 

89.29 

N  =  11 

13.77 

12.06 

100.58 

*4 

25.50 

31.75 

215.10 

N1  =  15 

32.52 

33.91 

36.71 

72.89 

N3  =  15 

N2  =  19 

6.1  Estimation  of  F,. 


Using  the  failure  times  given  in  Table  1,  we  obtain  via  (3.2)  an 
estimate  of  a  ,  a  =  16.9844  ;  note  that  this  estimate  is  different 
from  the  estimate  obtained  in  [9],  since  there,  all  the  data  of  [5]  were 
used.  We  then  rescale  the  data  in  Table  1  using  (3.3),  with  a  = 

16.9844  and  V^  =  20  ,  and  use  these  to  compute  the  empirical  distribu¬ 
tion  function  F^  ,  given  by  (3.4).  This  is  our  estimate  of  FQ  ,  and 


rr* 


12  - 


is  given  in  Figure  1.  The  median  of  these  rescaled  variables  is  easily 
verified  as  being  45572  minutes,  and  an  estimate  of  p  ,  obtained  from 
(3.5),  is  p  =  105998  minutes.  Other  s-percentiles  of  can  also  be 

A 

estimated  using  . 

6.2  Testing  Hypotheses 

A  plot  of  the  logarithms  of  the  60  rescaled  observations  men¬ 
tioned  in  Section  6.1  on  normal  probability  paper.  Figure  2,  suggests 
that  the  logarithms  of  the  failure  times  could  be  well  described  by  a 
normal  distribution. 

The  plot  of  Figure  2  suggests  that  we  test  for  normality  by 
taking  ,  i=l,2,3,4  ,  to  be  normal  distributions.  Our  estimates 

(u.,g.)  turn  out  to  be  (.902,  1.109),  (1.786,  1.525),  (2.228,  2.198), 

and  (3.821,  1.111);  D  is  computed  as  .4023.  From  the  table  of  criti¬ 
cal  values  given  in  Section  4.2  we  see  that  for  a  =  .05  ,  D*  =  .835  , 
and  since  D  <  D*  ,  we  have  no  reason  to  reject  the  null  hypothesis  that 
the  failure  times  of  Table  1  could  be  described  by  a  lognormal  distri¬ 
bution. 


The  approximate  95%  uniform  confidence  limits  for  F  are  ob¬ 
tained  as  Fq(0  t  ,835/v^O  ,  and  these  are  shown  by  the  dotted  lines  of 

Figure  1.  The  95%  confidence  interval  for  p  discussed  in  Section  5  is 
obtained  by  reading  the  ordinates  of  the  horizontal  line  at  .5  in  Figure 
1;  these  turn  out  to  be  32,320  and  64,000  minutes. 

Since  a  Weibull  distribution  for  times  to  failure  has  been  as¬ 
sumed  in  [5],  it  is  appropriate  that  we  test  for  this  distribution  as 
well.  We  thus  take  ,  i=l,...,4  to  be  an  extreme  value  distribu¬ 
tion  and  compute  (a^,3^)  as  (1.473,  1.154),  (2.531,  1.353),  (3.310, 
1.898),  and  (4.373,  .987);  we  obtain  D  as  .88.  From  the  table  of 
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Times  to  Failure  in 
Thousands  of  Minutes 


Figure  1.  Empirical  distribution  function  ?  of  the  60  rescaled 
observations  of  Table  1,  the  approximate  95%  uniform 
confidence  limits  for  F.,  and  a  confidence  interval 
for  V. 
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critical  values  given  in  Section  4.3  we  see  that  for  a  =  .05  ,  D*  = 
.87  .  Thus  D  is  significant  at  the  95%  level  of  significance,  but 
barely  so. 


Thus  it  appears  that  for  this  abstracted  set  of  data,  either  the 
lognormal  or  the  Weibull  distribution  would  provide  an  adequate  descrip¬ 
tion,  with  the  lognormal  having  a  slight  edge  over  the  Weibull.  In 
choosing  between  these  two,  one  should  remember  that  even  though  the 
formal  goodness  of  fit  tests  favor  the  lognormal,  the  Weibull  with  its 
monotone  failure  rate  may  be  a  more  appropriate  distribution  for  de¬ 
scribing  the  times  to  failure.  The  failure  rate  of  a  lognormal  distri¬ 
bution  first  increases  and  then  decreases,  and  this  may  not  be  realistic 
for  the  situation  considered  here. 


\ 


'*.4 
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